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Cn ' Abstract. Let K be the number field determined by a monic irreducible 

polynomial f{x) with integer coefficients. In previous papers we parameter- 
^ ^ ized the prime ideals of K in terms of certain invariants attached to Newton 

polygons of higher order of the defining equation f{x). In this paper we show 

\K\ ' how to carry out the basic operations on fractional ideals of K in terms of 

these constructive representations of the prime ideals. From a computational 
perspective, these results facilitate the manipulation of fractional ideals of K 
avoiding two heavy tasks: the construction of the maximal order of K and the 

tn ' factorization of the discriminant of f{x). The main computational ingredient 

y-^ ' is Montes algorithm, which is an extremely fast procedure to construct the 

^^ prime ideals. 

-)— > 

Introduction 

-Y-\ ! Let if be a number field of degree n and Z^ its ring of integers. An essential 

^ ' task in Algorithmic Number Theory is to construct the prime ideals of K in terms 

\^ \ of a defining equation of K , usually given by a monic and irreducible polynomial 

'^ ' f{x) e T^x]. The standard approach to do this, followed by most of the algebraic 

manipulators like Kant, Pari, Magma or Sage, is based on the previous computation 
of an integral basis of Z^- This approach has a drawback: one needs to factorize 

1^ \ the discriminant, disc(/), of /(x), which can be a heavy task, even in number fields 

^■"■^ ■ of low degree, if f{x) has large coefficients. 

In this paper we present a direct construction of the prime ideals, that avoids 
the computation of the maximal order of K and the factorization of disc(/). The 
following tasks concerning fractional ideals can be carried out using this construc- 
tion: 

(1) Compute the p-adic valuation, Vp : K* -^ Z, for any prime ideal p of K. 

(2) Obtain the prime ideal decomposition of a fractional ideal. 

(3) Compute a two-element representation of a fractional ideal. 

(4) Add, multiply and intersect fractional ideals. 

(5) Compute the reduction maps, Jjk -^ Z/^/p. 

(6) Solve Chinese remainders problems. 
Moreover, along the construction of a prime ideal p, lying over a prime number 

p, a Zp-basis of the ring of integers of the local field Kp is obtained as a by-product. 
Hence, from the prime ideal decomposition of the ideal pLk we are also able to 
derive the resolution of another task: 
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(7) Compute a p-integral basis of K. 

For a given prime number p, the prime ideals of K lying above p are in one- 
to-one correspondence with the irreducible factors of f{x) in Zp[a;] |Hen08j . In 
}HN08I we proved a series of recurrent generalizations of Hensel lemma, leading to a 
constructive procedure to obtain a family oi f- complete types, that parameterize the 
irreducible factors of f{x) in Zp[x]. A type is an object that gathers combinatorial 
and arithmetic data attached to Newton polygons of f{x) of higher order, and 
an /-complete type contains enough information to single out a p-adic irreducible 
factor of f{x). In |GMN08] we described Montes algorithm, which optimizes the 
construction of the /-complete types; it outputs a list of /-complete and optimal 
types that parameterize the prime ideals of K lying above p, and contain valuable 
arithmetic information on each prime ideal. All these results were based on the 
PhD thesis of the second author |Mon99| . The algorithm is extremely fast in 
practice; its complexity has been recently estimated to be 0{n^^'^6 + n^^'^d^^'^), 
where .5 = log(di sc(/)) [FvTO] . 

In |GMN09b| we reinterpreted the invariants stored by the types in terms of 
the Okutsu polynomials attached to the p-adic irreducible factors of f{x) |Oku82) . 
Suppose t is the /-complete and optimal type attached to a prime ideal p, cor- 
responding to a monic irreducible factor fp{x) G Zp[x]; then, the arithmetic in- 
formation stored in t is synthesized by two invariants of /p (x) : an Okutsu frame 
[(j)i(x), . . . ,(j)r{x)] and a Montes approxim,ation (l)p{x) (cf. loc.cit.). The monic 
polynomials (/)i, . . . , 0^, </>p have integer coefficients and they are all irreducible over 
Zp [x] ; the polynomial 4>p (x) is "sufficiently close" to /p {x) . We say that 

P = [p;<Pi,---,(t>r,<Pp], 

is the Okutsu-Montes representation of the prime ideal p. Thus, from the compu- 
tational point of view, p is structured in r -|- 1 levels and at each level one needs 
to compute (and store) several Okutsu invariants that are omitted in this nota- 
tion. This computational representation of p is essentially canonical: the Okutsu 
invariants of p, distributed along the different levels, depend only on the defining 
equation f{x). These invariants provide a rich and exhaustive source of informa- 
tion about the arithmetic properties of p, which is crucial in the computational 
treatment of p. 

From a historical perspective, the sake for a constructive representation of ideals 
goes back to the very foundation of algebraic number theory. Kummer had the 
insight that the prime numbers factorize in number fields into the product of prime 
"ideal numbers" , and he tried to construct them as symbols [p ; </»] , where (^(x) is a 
monic lift to Z[a;] of an irreducible factor of f{x) modulo p. Dedekind showed that 
these ideas led to a coherent theory only in the case that p does not divide the index 
?(/) := [Ijk'- "^Wi/ifix))). This constructive approach could not be universally 
used because there are number fields in which p divides the index of all defining 
equations jDed78| . Fortunately, this obstacle led Dedekind to invent ideal theory as 
the only way to perform a decent arithmetic in number fields. Ore, in his Phd thesis 
J0re23j . tried to regain the constructive approach to ideal theory. He generalized 
and improved the classical tool of Newton polygons and showed that under the 
assumption that the defining equation is p-regular (a much weaker condition than 
Dedekind's condition p \ i{f)), the prime ideals dividing p can be parameterized 
as p = [p ; </), (j)p] (in our notation), where 4>p{x) G Z[x] is certain polynomial whose 
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0-Newton polygon is one-sided and the residual polynomial attached to this side is 
irreducible (cf. section [1]). The contribution of |Mon99| was to extend Ore's ideas 
in order to obtain a similar construction of the prime ideals in the general case. 

The aim of this paper is to show how to use this constructive representation of 
the prime ideals to carry out the above mentioned tasks (l)-(6) on fractional ideals 
and to compute p-integral bases. The outline of the paper is as follows. In section 
[T]we recall the structure of types, we describe their invariants, and we review the 
process of construction of the Okutsu-Montes representations of the prime ideals. 
In section [2] we show how to compute the p-adic valuation of K with respect to a 
prime ideal p; this is the key ingredient to obtain the factorization of a fractional 
ideal as a product of prime ideals (with integer exponents). The operations of 
sum, multiplication and intersection of fractional ideals are trivially based on these 
tasks. In section[3]we show how to find integral elements ap € %k such that p is the 
ideal of 'Lk generated by p and ap ; this leads to the computation of a two-element 
representation of any fractional ideal. In section|4l we show how to compute residue 
classes modulo prime ideals and we design a Chinese remainder theorem routine. 
Section [5] is devoted to the construction of a p-integral basis. 

We have implemented a package in Magma that performs all the above mentioned 
tasks; in section IH] we present several examples showing the excellent performance 
of the package in cases that the standard packages cannot deal with. Our routines 
work extremely fast as long as we deal only with fractional ideals whose norm may 
be factorized. Even in cases where disc(/) may be factorized and an integral basis 
of Z/f is available, our methods work faster than the standard ones if the degree of 
K is not too small. Mainly, this is due to the fact that we avoid the use of linear 
algebra routines (computation of Z-bases of ideals, Hermite and Smith normal 
forms of n X n matrices, etc.), that dominate the complexity when the degree n 
grows. Finally, in section[7]we make some comments on the apparent limits of these 
Montes' techniques: they are not yet able to test if a fractional ideal is principal. 
We also briefiy mention how to extend the results of this paper to the function field 
case and the similar challenges that arise in this geometric context. 

Notations. Throughout the paper we fix a monic irreducible polynomial ](x) G 
1\x\ of degree n, and a root G Q of /(x). We let K ~ Q{0) be the number field 
generated by 6, and Zk its ring of integers. 

1. Okutsu-Montes representations of prime ideals 

Let p be a prime number. In this section we recall Montes algorithm and we 
describe the structure of the /-complete and optimal types that parameterize the 
prime ideals of K lying over p. The results are mainly extracted from [HN08^ (HN 
standing for "Higher Newton") and |GMN08| . 

Given a field F and two polynomials ip{y), 4'{y) G F[y], we write fijj) '^ ''Piu) 
to indicate that there exists a constant c ^ F* such that ip{y) = cip{y). 

1.1. Types and their invariants. Let v. Qp — > Q be the canonical extension of 
the p-adic valuation of Qp to a fixed algebraic closure. We extend v to the discrete 
valuation vi on the field Qp{x), determined by: 

vi : Qp[x] — >ZLI {oo}, vi{bo + ■ ■ ■ + Wx'') := mm{v{bj), < j < r}. 
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Denote by Fq := GF{p) the prime field of cfiaracteristic p, and consider the 0-th 
residual polynomial operator 



i?o: Zp[x] -^ ¥o[y], g{x) ^ g{y)/p-^(a), 

where, : 1'p[y] — >■ Fo[y], is the natural reduction map. A type of order zero, t = 
ipoiy)i is just a monic irreducible polynomial ipa{y) € ]Fo[y]- A representative oft is 
any monic polynomial (/'i(x) € Z[a;] such that i?o('/'i) = ipo- The pair {(j)i,vi) can be 
used to attach a Newton polygon to any nonzero polynomial g{x) G Qp[x] . If g{x) — 
^j>g as(a;)0i(a;)* is the (/)i-adic development of g{x), then Ni{g) := N^j^^y^{g) 
is the lower convex envelope of the set of points of the plane with coordinates 
{s,vi{as{x)(t>i{xy)) [HN081 Seel]. 

Let Ai G Q~ be a negative rational number, Ai = —hi/ei, with hi,ei po- 
sitive coprime integers. The triple (0i, wi, Ai) determines a discrete valuation V2 on 
Qp{x), constructed as follows: for any nonzero polynomial g{x) £ Zp[x], take a line 
of slope Ai far below Ni{g) and let it shift upwards till it touches the polygon for 
the first time; if H is the ordinate at the origin of this line, then V2{g{x)) = eiH, by 
definition. Also, the triple (0i,wi, Ai) determines a residual polynomial operator 

i?i := R^,,v,,Xi ■■ ^pN ^ ]Fi[y], Fi := ¥o[y]/{i;oiy)), 
which is a kind of reduction of first order of g{x) HN08[ Def.1.9]. 

Let V'i(y) G '^i[y] be a monic irreducible polynomial, ipi{y) 7^ y. The triple 
t = (01 (x) ; Ai, 7/11(2/)) is called a type of order one. Given any such type, one can 
compute a representative of t; that is, a monic polynomial 4'2ix) € Z[x] of degree 
ei degf/"! deg^i, satisfying Ri{(j)2){y) ~ '^i{y)- Now we may start over with the 
pair (02, ^^2) and repeat all constructions in order two. 

The iteration of this procedure leads to the concept of type of order r |HN08[ 
Sec. 2]. A type of order r > 1 is a chain: 

t = [cj)i{xy,\i,cj)2{x)]-- ■ ;Ar_i,0r(x);Ar,V'r(y)), 
where <f)i{x), . . . , 4>r{x) are monic polynomials in Z[a:] that are irreducible in 'Lp[x\, 
Ai, . . . , Ar are negative rational numbers, and "iprijj) is a polynomial over certain 
finite field F^ (to be specified below), that satisfy the following recursive properties: 

(1) ipiix) is irreducible modulo p. We define V'o(2/) := Ro{4'i){y) G ^o[y], 

¥,^¥o[y]/{My))- 

(2) For all 1 < z < r, A'^i(0i+i) := N^.^y.{(j)i+i) is one-sided of slope Ai, and 

Ri{4>t+i){y) ■■= i?0,,i,^.A. (0i+i)(j/) ~ ipi{y), for some monic irreducible poly- 
nomial V»(2/) G Hy]- We define F,+i = ¥,[y]/{,p,{y)). 

(3) Tpriy) S ^r[y] is a monic irreducible polynomial, ipriy) 7^ y- 

Thus, a type of order r is an object structured in r levels. In the computational 
representation of a type, several invariants are stored at each level, 1 < i < r. The 
most important ones are: 

4>i{x), monic polynomial in Z[a;], irreducible over Zp[x] 

rui, deg(j)i{x), 

Vi{(j)i) non- negative integer, 

Ai = —hi/ei, hi,ei positive coprime integers, 

(!-{,(-[, a pair of integers satisfying dihi — i'^ei = 1, 

ipiiy), monic irreducible polynomial in Fi[y], 

fi degtp.iy), 

Zi the class of y in Fi+i, so that ipi{zi) = 0. 
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Take /o := degi/'o, and let za G Fi be the class of y, so that iJjq{zo) = 0. Note 
that rrii = (/o/i • • • /j-i)(ei • • • e^-i), Fj+i = F^^i], and diniFo Fj+i = /o/i • • ■ /». 
The discrete valuations ui, . . . , u,-|-i on the field Qp{x) are essential invariants of 
the type. 

Definition 1.1. Let g{x) G Zp[a;] be a monic separable polynomial, and t a type of 
order r > 1. 

(1) We say that t divides g{x) (and we write 1 1 g{x)), if tpriv) divides Rj.[g){y) 
in¥r[y]. 

(2) We say that t is g-complete if OTd^,^,{Rr{g)) = 1. In this case, t singles 
out a monic irreducible factor gt{x) G Zp[x] of g{x), uniquely determined by the 
property Rr{gt){y) ~ ipriv)- If Kt *s the extension of Qp determined by gtix), then 

eiKt/Qp) = ei ■ • ■ e,, fiKt/Qp) - /o/i ■ ■ ■ fr- 

(3) A representative of t is a monic polynomial (l)r+i{x) G Z[a;], of degree 
rrir+i = erfr'TTir such that Rr{(t>r+i)iy) ~ i^riv)- This polynomial is necessarily 
irreducible in Zp[x]. By the definition of a type, each 0i+i(a:) is a representative of 
the truncated type of order i 

Trunci(t) := (0i(a;); Ai, ^2(2;); ■ • • ;Xi^i,(j>i{x);Xi,ipi{y)). 

(4) We say that t is optimal if mi < • • • < rur, or equivalently, if Cifi > 1, for 
all I < i < r. 

Lemma 1.2. Let t be a type of order r. Then, Vj{4'i) = {mi/mj)vj{(f)j), for all 
j < i < r. 

Proof. Let (t)i{x) = J2s>o'^s{x)(j)j{x)^ be the (j)j-adic development oi(j)i. By |HN08| 
Lem.2.17], Vj{(j)i) = mi\is>Q{vj[as(f)^)}- Now, Nj{(j)i) is one-sided of slope A^-, 
because (pi is a polynomial of type TrunCj(t) [HN08[ Def.2.1-|-Lem.2.4]. Since 
the principal term of the development is <^™''™j^ -y^re get Vj{4)i) = Uj(0™ ) = 
{mi/mj)vj{(j)j). a 

1.2. Certain rational functions. Let t be a type of order r. We attach to t seve- 
ral rational functions in Q{x) |HN08| Sec. 2. 4]. Note that Vi{ipi) is always divisible 
by ei_i IHNOBI Thm.2.11]. 

Definition 1.3. Let 7ro(a;) = 1, 7ri(x) —p. We define recursively for all 1 < i < r: 



These rational functions can be written as a product of powers of p, (j)i{x), . . . , (j)r{x), 
with integer exponents. 

Notation. Let '^{x) = p"°<j>i{x)'^'^ ■ ■ ■ (psix)""" G Q{x) be a rational function which 
is a product of powers of p, (pi, ... ,(j)s, with integer exponents. We denote: 

log* = (no,...,n,)GZ-'+i. 

The next result is inspired in |HN08[ Cor. 4. 26]. 

Lemma 1.4. Let F{x) G Zp[a;] be a monic irreducible polynomial divisible byt, and 
let a eQp be a root ofF{x). For some 1 < s <r, let'^{x) = p"°(j)i{x)"^ ■ ■ -(psixY^" 
be a rational function in Q,{x), such that u(^(a)) = 0. Then, 

*(a;) =71(2;)*' ■■■IsixY", 
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for certain integer exponents ti, . . . ,ts G Z, which can be computed by the following 
recursive procedure: 

vector=(no, . . . , rig) 
for i=s to 1 by -1 do 
ti=vector [i] /e^ 
vector=vector-ti log ji 
end for 

Proof. By (HN081 (17)] and |HN08l Cor.3.2]: 

log$s = ( , 1), log^, = ( ,0), log7, = ( , e,) e Z'+\ 

w(0s(a)) ^y^eifi---es^ifs~i — , v{js{a)) =0. 



i=l 



ei • ■ • Cj 



Since i'(\l/(a)) = 0, the formula for v((j)s{a)) shows that esjris- Thus, we can replace 
\l/(a;) by \l/(a;)7s '^''^' and iterate the argument. Since u(7s(a)) = 0, the new ^(x) 
satisfies w(^(a)) = as well, and the s-th coordinate of log$ is zero. At the 
last step (s = 1), we get 'i'{x) — p"o0"\ with u'q + n[{hi/ei) = 0. Then, clearly 
*(a;) =7i(x)"'i/'=i. □ 

1.3. Montes algorithm and the secondary invariants. Let f{x) G Z[x] be a 
monic irreducible polynomial. At the input of the pair {f{x),p), Montes algorithm 
computes a family ti, . . . ,ts of /-complete and optimal types in one-to-one cor- 
respondence with the irreducible factors ft^{x), . . . , ftS^) of f{x) in "Lplx]. This 
one-to-one correspondence is determined by: 

(1) For all 1 < i < s, the type t,j is /t.-complete. 

(2) For all j ^ «, the type tj does not divide fti{x). 

The algorithm starts by computing the order zero types determined by the irre- 
ducible factors of f{x) modulo p, and then proceeds to enlarge them in a convenient 
way till the whole list of /-complete optimal types is obtained |GMN08) . 

With regard to the computation of generators of the prime ideals and Chinese 
remainder multipliers, the algorithm is slightly modified to compute and store some 
other (secondary) invariants at each level of all types t considered by the algorithm: 

RefinementSi, a list of pairs [(j){x), A], where </> is a representative of 

TrunCi_i(t) and A a negative slope, 

Ui, a nonnegative integer called the height, 

Quotj, a list of ei polynomials in 1i[x], 

log<l>i, a vector {uq, . . . ,ni) € Z'+^, 

logTTi, a vector (no, . . . ,rii_i,0) G Z''+^, 

log7i, a vector {no, . . . ,ni) 6 Z*+^. 

Let us briefiy explain the flow of the algorithm and the computation of these 
invariants. Suppose a type of order i — 1 dividing f{x) is considered, 

t = {4)i{x);\i,(l)2{x);--- ; Ai_2,0j-i; Ai_i,V'i-i(y))- 

A representative 0i(x) is constructed. Suppose that either i = 1 or mi < ■ ■ ■ < mi. 
Let I = ord^._i i?i_i(/). If t is not /-complete (£ > 1), it may ramify to produce 
new types, that will be germs of distinct /-complete types. To carry out this 
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ramification process we compute simultaneously the first £ + 1 coefficients of the 
(/)i-adic development of f{x) and the corresponding quotients: 

f{x) = (p.i{x)qi{x)+ao{x), 
Qx qi{x) = (j)i{x)q2{x) +ai{x), 

qe{x) = (l)i{x)qi+i{x) +ae{x). 

The Newton polygon of i-th order of f{x), Ni{f), is the lower convex envelope 
of the set of points {s,Vi{as4'f)) of the plane, for all s > 0. However, we need 
to build up only the principal part of this polygon, N^^ (f) — N7 ^. (/), formed 
by the sides of negative slope of Ni{f). By [HN08, Lem.2.17], this latter polygon 
is the lower convex envelope of the set of points {s,Vi{as(j)f)) of the plane, for 
< s < £. For each side of slope (say) A of N~{f), the residual polynomial 
R^{f){y) ~ R'Pi,vi,\{f){y) & ^i[y] is computed and factorized into a product of 
irreducible factors. The type t branches in principle into as many types as pairs 
(A, ■0(2/)), where A runs on the negative slopes of N~{f) and tp{y) runs on the 
different irreducible factors of Rx{f){y)- If one of these branches 

tA,^ := (01 (a;); Xi,(j)2{x); • • • ; Aj_i, 0,(x); A, il^iy)), 

is /-complete, we store this type in an specific list and we go on with the analysis 
of other branches. Otherwise, we compute a representative (t)x,ip{x) of t^^^. Let e 
be the least non-negative denominator of A and / — deg ip. Then we proceed in a 
different way according to e/ = 1 or e/ > 1. 

If ef > 1, then deg 0a,?/; > "t-i, so that (pi+i :— 0a, V' ™-&y be used to enlarge t^^^ 
into several optimal types of order i -I- 1. We store the invariants 

0i, rui = deg0i, Vi{<pi), Xi = A, hi, e^ = e, ii, i\, ipi ^ip, fi^ /, Zi, 
Ui, Quotj, log$i, logTTi, log7i 

at the j-th level of t^,^), and then we proceed to enlarge the type. The invariant 
Vi{(j)i) is recursively computed by using |HN08[ Prop.2.7+Thm.2.11]: 

v(6) = I '^' if i = 1, 

'■^^'' \ ei_i/i_i(ei_iWi_i(0i_i) -|-/ii_i), ifi>l. 

Let Si be the abscissa of the right end point of the side of N,^' (/) of slope A; the 
secondary invariants are computed as: 

Ui = Vi{asi), 

Quot, = [1, gs._i(x), . . . , qs,-e+i{x)], 
(2) log*! = (no,...,ni_i,l), where (no, . . . ,nj_i) = -(wi(0i)/ei_i) logTr^.i, 
logTTi =£i_ilog$i_i -£-_ilog7ri_i, 
log 7i = Ci log $i - hi log TTi . 

For i = 1 we take log$i = (0, 1) and logTTi = (1, 0). Note that all these secondary 
invariants depend only on A and not on ip. They are computed only once for each 
side of N~ (/) and then stored in the different optimal branches that share the same 
slope. The type t;^^^ of order i is then ready for further analysis. 

If e/ = 1 then deg 4>\^^ = rrii, so that the enlargements of t^^^ that would result 
from building up an {i + l)-th level from 0>^^ would not be optimal. In this case. 
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we replace t\^^ by all types 

obtained by enlarging t with i-th levels deduced from the consideration of 0a.i/j as a 
new (and better) representative oft: 0,^(x) :— 4>\^^{x), but taking into account only 
the slopes A' of N^'.,Vi (/) satisfying A' < A; this is called a refinement step [GMNOSi 
Sect. 3. 2]. If the total number of pairs (A, ip) is greater than one, we append to the 
list Refinements^ of all types t^, ^, the pair [(f)i(x), A]. Some of these new branches 
*A' lb' '-'^ order i may be /-complete, some may lead to optimal enlargements and 
some may lead to further refinement at the i-th level. Thus, in general, the list 
RefinementSi of a type t is an ordered sequence of pairs: 

Refinements, = [[cj)^\ X^\ • • • , [cl^^\xh], 
reflecting the fact that along the construction of t, s or more successive refine- 
ment steps occurred at the z-th level. All polynomials (j)\ are representatives of 
TrunCi_i(t), and the slopes A^ grow strictly in absolute size: |A^ | < • ■ • < |A^* |. 
We recall that we store only the pairs [4>l , X\ ] corresponding to a refinement step 
that occurred simultaneously with some branching. For instance, if Af"^) (/) has 

only one side with integer slope A- G Z (i.e. e = 1) and the corresponding residual 
polynomial is the power of an irreducible polynomial of degree one (/ = 1), then 
the pair [4>i \Xl ] is not included in the list Refinements^. 

After a finite number of branching, enlargement and/or refinement steps, all 
types become /-complete and optimal. If t is an /-complete and optimal type of 
order r, then the Okutsu depth of /t(a;) is |GMN09b|. Thm.4.2]: 

(3) R = 

The invariants fi+i, hi, Ci, fi at each level I < i < R are canonical (depend 
only on f{x)) |GMN09bl Cor. 3. 7]. On the other hand, the polynomials 4>i{x), ipiiv) 
depend on several choices, some of them caused by the lifting of elements of a 
finite field to rings of characteristic zero. However, the sequence [(pi, . . . , 4'ji\ is an 
Okutsu frame of ft{x) [GMN09bl Sec. 2]; in the original terminology of Okutsu, the 
polynomials (/)i, . . . , (/i/j are primitive divisor polynomials of ft{x) [Oku82| . 

1.4. Montes approximations to the irreducible p-adic factors. Once an /- 
complete and optimal type t of order r is computed, Montes algorithm attaches to 
it an (r + l)-level that carries only the invariants: 

/ .s 4'r+l, fnr+1, Vr+l{4>r+l), A,.+ i = — /ir+l, 6^+1 = 1, V-'r+l- fr+1 = I7 

Zr+1, l0g$r+l, logTTr+l, log7r+l. 

The polynomial (jjr+iix) is a representative of t. The invariants A^+i, 4'r+i are 
deduced from the computation of the principal part of the Newton polygon of 
(r -I- l)-th order of f{x) (which is a single side of length one) and the corresponding 
residual polynomial Rr+i{f){y) G Fr+i[y] (which has degree one). 
If p is the prime ideal corresponding to t, we denote 

fp{x) := ft{x) e Zp[x], 0p(x) := 0,.+i(x) £ Z[x], Fp := F,+i, 
tp := ((/)i;Ai,^2;--- ,(/)r;Ar,^p;A,.+i,i/'r+i), 



r, 




if Crfr > 1, 


r - 


- 1, 


if Crfr - 1. 
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and we say that p = [p;(f)i, . . . ,(l)r, (f>p] is the Okutsu-Montes representation of p. 

Note that tp is still an /-complete type of order r + 1, but it may eventually be 
non-optimal because m^+i = m^, if Crfr = 1- 

The polynomial 4>p {x) is a Montes approximation to the p-adic irreducible factor 
fp{x) [GMNO9R Sec. 4.1]. Several arithmetic tasks involving prime ideals (tasks (1), 
(3), (5), (6) and (7) from the list given in the Introduction) require the computation 
of a Montes approximation with a sufficiently large value of the last slope | A^+i | = 
hr+i- This can be achieved by applying a finite number of refinement steps at the 
(r + l)-th level, to the type tp of order r -t- 1, as described in |GMN09bl Sec. 4. 3]. 
This procedure has a linear convergence. In |GNP10j a more efficient single-factor 
lift algorithm is developed, which is able to improve the Montes approximations to 
fp{x) with quadratic convergence. 

2. p-ADIC VALUATION AND FACTORIZATION 

In section 12.11 we compute the p-adic valuation, Vp : K* -^ Z, determined by a 
prime ideal p, in terms of the data contained in the Okutsu-Montes representation 
of p. In section [221 we describe a procedure to find the prime ideal decomposition 
of any fractional ideal of K. From the computational point of view, this procedure 
is based on three ingredients: 

(1) The factorization of integers. 

(2) Montes algorithm to find the prime ideal decomposition of a prime number. 

(3) The computation of Vp for some prime ideals p. 

The routines (2) and (3) run extremely fast in practice (see section [5]). 

It is well-known how to add, multiply and intersect fractional ideals once their 
prime ideal factorization is available. We omit the description of the routines that 
carry out these tasks. 

From now on, to any prime ideal p oi K we attach the data tp, fp{x), (f>p{x), Fp, 
described in section FOl Also, we choose a root 9p G Qp of fp{x), we consider the 
local field Kp = Qp{Op), and we denote by Zk^ the ring of integers of A'p. 

2.1. Computation of the p-adic valuation. Let p be a prime number and 
v: Q — > Q, the canonical p-adic valuation. Let p be a prime ideal of K lying 
above p, corresponding to an /-complete type tp with an added (r -I- l)-th level, as 
indicated in section [1^41 We shall freely use all invariants of tp described in section 
[H By item 2 of Definition 1 1.1 1 we know that 

e(p/p) = ei • ■ • e^, f{p/p) = /o/i ■ ■ ■ fr- 

The residue field Zkp /p'^Kp can be identified to the finite field 

Fp := ¥r+i = Fo[zo,2;i, . . . ,Zr]. 

More precisely, in |HN08 '. (27)] we construct an explicit isomorphism 



(5) 7:Fp ^^ ZkJpZkp, zq <-^ 0p, zi ^ jiiOp), . . . , Zr ^ jriOp), 

where we indicate by a bar the canonical reduction map, Zk^ — ;■ Zk^ /pZ^'p . We 
denote by Iredp : Zk^ — > Fp , the reduction map obtained by composition of the 
canonical reduction map with the inverse of the isomorphism ([5]). 

(6) Iredp : Zk^ -^ Z^, /pZ^, ^ Fp. 
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Consider the topological embedding Lp: K ^^ Kp determined by sending 6 to 9p. 
We have: Vp{a) — e(p/p)v{ip{a)), for all a G K. In particular, for any polynomial 
g{x) e Z[x], 

(7) vp{gie))^e{p/p)v{g{ep)). 

Any a € K* can be expressed as a = {a/b)g{9), for some coprime positive 
integers a, b and some primitive polynomial g{x) € Z[a;]. By (O, 

vp{a)^e{p/p){v{g{ep))+v{a/b)). 

Thus, it is sufficient to learn to compute v{g{6p)). The condition v{g{9p)) = is 
easy to check |GMN09bl Lem.2.2]: 



(8) 



vig{0p)) = 



V^ot-RoCff)- 



If V'o I Ro (g) J the computation oiv{g{9p)) can be based on the following proposition, 
which is easily deduced from |HN081 Prop.3.5] and jHN08[ Cor.3.2]. 

Proposition 2.1. Let p, fp{x), 0p be as above. Let t be a type of order R dividing 
fp{x), and let g{x) £ Z[a;] be a nonzero polynomial. For any I < i < R, take a line 
Lxi of slope Xi far below Ni{g), and let it shift upwards till it touches the polygon 
for the first time. Let S be the intersection of this line with Ni{g), let {s,u) be the 
coordinates of the left end point of S, and let H = u + s\Xi\ be the ordinate at the 
origin of this line. Then, 

(1) v{g{6p)) > H/ei ■ ■ -e-i-i, and equality holds if and only i/TrunCi(t) \ g(x). 

(2) If equality holds, then v{g(6p)) = v{'^i{dpY'ni(6pY) and 



Iredn 



9{0p 



^^{epY^^{epY 



R^{g){z^)^0. 



D 



Figure 1 shows that the segment S may eventually be reduced to a point. In 
this case, the residual polynomial Ri{g){y) is a constant |HN08[ Def.2.21], so that 
TrunCi(t) \ g{x) automatically holds. 





Figure 1 

We may compute Vp{g{9)) ~ e{p/p)v{g{9p)) by applying Proposition 12.11 to the 
type tp. If for some 1 < i < r + 1, the truncated type TrunCi(tp) does not divide 
g{x), we compute v{g{9p)) as indicated in item 1 of this proposition. Nevertheless, 
it may occur that TrunCi(tp) divides g{x) for all 1 < i < r + 1 (for instance, if 
g{x) is a multiple of <^p(a;) — 0r+i(a;)). In this case, we compute an improvement of 
the Montes approximation 0p (x) by applying the single- factor lift routine }GNP10) ; 
then, we replace the (r + l)-th level of tp by the invariants Q determined by the 
new choice oi(j)r+iix) — (fipix), and we test again if tp = TrunCr+i(tp) divides g{x). 
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If tp divides g{x), then 4>pix) is simultaneously close to ap-adic irreducible factor 
of f{x) and to a p-adic irreducible factor of g(x); hence, if f{x) and g{x) do not 
have a common p-adic irreducible factor, after a finite number of steps the renewed 
type tp will not divide g{x). On the other hand, if f{x) and g{x) have a common 
p-adic irreducible factor, they must have a common irreducible factor in Z[x] too; 
since f{x) is irreducible, necessarily f{x) divides g{x) and g{9) — 0. 

We may summarize the routine to compute Vp{a) as follows. 

Input: a G K* and a prime ideal p determined by a type tp of order r + 1. 
Output: Vp{a). 

1. Write a = f5(^), with a, 6 coprime integers and g{x) G Z[a:] primitive. 

2. Compute i^ — v{a/h). 

3. if -i/jo t Ro{g) then return Up(a) = e(p/p)z^. 

4. for i = 1 to r + 1 do 

compute Nf{g), Ri{g), and the ordinate H of Proposition [Ol 
if ■01 1 -Ri(5) then return Wp(a) = e{p/p){{H/ei • ■ ■ e^-i) + v). 
end for. 

5. while V'r+i I ^r+i(5) do 

improve (f)p and compute the new values A^+i, i^r+i- 
compute N^_^-^^{g), Rr+i{g), and the ordinate H of Proposition [^TTJ 
end while. 

6. return Vp{a) ~ e{p/p){{H/ei ■ ■ ■ e,-) + i^). 

2.2. Factorization of fractional ideals. For any a G K*, the factorization of 
the principal ideal generated by a is 

p 

Let a = {a/b)g{6), for some positive coprime integers a,b and some primitive 
polynomial g{x) G Z[x]. Then, Vp{a) = for all prime ideals p whose underlying 
prime number p does not divide the product abNx/q{g{9)) = a6Resultant(/, g). 
Also, if p is a prime ideal of K and o, b are fractional ideals, we have 

Vp{a+ b) = min{t;p(a),i;p(b)}. 

Thus, the p-adic valuation of the fractional ideal a generated by ai, . . . , a,„ G K* 
is: tip (a) = m.mi<i<rn{vp{ai)}. 

After these considerations, it is straightforward to deduce a factorization routine 
of fractional ideals from the routines computing prime ideal decompositions of prime 
numbers and p-valuations of elements of K* with respect to prime ideals p. 

Input: a family ai, . . . , am S K* of generators of a fractional ideal a. 
Output: the prime ideal decomposition a = Yip P'^" ■ 

1. For each I < i < m, write at = {ai/bi)gi{9), with ai,bi coprime integers and 
gi{x) G Z[a;] primitive; then compute Ni = N^/q (17^(6')). 

2. Compute N = gcd{aiNi, . . . , amNm) and M = lcm(6i, . . . , bm)- 

3. Factorize N and M and store all their prime factors in a list V. 

4. For each p G P apply Montes algorithm to obtain the prime ideal decomposition 
of p, and for each p|p, take Op = mini<i<,„{t;p(Q;i)}. 

5. Return the list of pairs [p, Op] for all p with ap ^ 0. 
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The bottleneck of this routine is step 3. We get a fast factorization routine in 
the number field K, as long as the integers iV, M attached to the ideal o may be 
easily factorized. 

3. Computation of generators 

In |GMN08l Sec. 4] we gave an algorithm to compute generators of the prime 
ideals as certain rational functions of the (/)-polynomials. Some inversions in K, 
one for each prime ideal, were needed. These inversions dominated the complexity 
of the algorithm, and they were a bottleneck that prevented the computation of 
generators for number fields of large degree. 

In sections [3.11 and [3^ we construct a two-element representation of prime ideals, 
which does not need any inversion m. K . As a consequence, this construction works 
extremely fast in practice even for number fields of large degree (see section [5]). In 
section [3751 we easily derive two-element representations of fractional ideals. 

For any prime ideal p oi K we keep the notations for tp, fp{x), <^p(a;), Fp, 0p, 
Kp , Zxp , as introduced in section [51 

3.1. Local generators of the prime ideals. 

Definition 3.1. A pseudo-generator of a prime ideal p of K is an integral element 
TT G liK such that Wp(7r) = 1. 

Let p be a prime number, and let p — [p] 0i, . . . , (pr, (j)p\ be a prime ideal factor 
of pZ/j-, corresponding to an /-complete type tp with an added (r -I- l)-th level, as 
indicated in section [L4l In this section we show how to compute a pseudogenerator 
of p from the secondary invariants m^, Quotj, of tp, for 1 < i < r, computed along 
the flow of Montes algorithm as indicated in ([2]) . 

For each level 1 < i < r, let us denote: 

QuOtj = [Qifi[x), . . . ,Qi^e,~l{x)]. 

Recall that Qi,o{x) — 1, and for < j < e^, the polynomial Qi,j{x) G Z[x] is the 
{si — j)-th quotient of the (j)i-a,dic development of f{x) (cf. dl])), where Si is the 
abscissa of the right end point of the side of slope A^ of N~{f). Also, let us define 

H- n H- u,+j{\X,\+Vi{(l},)) wp,/-/ 

Hi^a = 0, Hij = , VO < J < Ci. 

ei ■■■Bi-i 

Proposition 3.2. For each level I < i < r and suhindex < j < ei: 

(1) Vcf{Qi j{9)) > e{q/p)Hi j, for all prime ideals q \ p. 

(2) Vp{Ql,{9))=eip/p)H-^,. 

Proof. Item 1 being proved in |GMN09a[ Prop. 10], let us prove item 2. Fix a 
level 1 < i < r and a subindex < j < Ci. Let ii = ord^;_^ i?i_i(/), and let 
f{x) = J2s>o "s'^i ^® ^^^ (j)i-a,dic development of f{x). The Newton polygon of i-th 
order of f{x), Ni{f), is the lower convex envelope of the cloud of points (s, Wi(as0f )), 
for aU s > 0. The principal part N~{f) is equal to Ni{f) n ([0, £i\ x R); the typical 
shape of this polygon is illustrated in Figure 2. Let S\. be the side of slope Xi of 
this polygon, and Si the abscissa of the right end point of S\- . 
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Figure 2 

For any < t < i?i, let qt{x) be the t-ih quotient of the 0i-adic development (see 
(P)). We have f{x) = qt{x)(f)i{xy + rt{x), with 



rt{x) 



E 

0<s<t 



as{x)(l)t{xy , qt{x)(f)i{xY 



'^as{x)(j)i{xy 

t<s 



Hence, if to is the smallest abscissa of a vertex of Ni{f), such that ip > t, we have 

N,{qt(f>l) n ([io, w) X R) = N,{f) n {[to, w) X K) . 

Recall that Qij{x) — qt{x), for t — Si — j; for this value of t we have to — Si (see 
Figure 2). On the other hand, all points in the cloud (s, Ui(as(/>f)), for s^—Ci < s < Si 
lie strictly above S\- , because the point on 5*^^ with integer coordinates and closest 
to the right end point has abscissa Si — e^. Hence, the line of slope A^ that first 
touches Ni{qs.^j(j)^'~-') from below is the line containg S\-. This line has ordinate 
at the origin (see Figure 2): 

H ^ Vi {asi4i^y) + Si\Xi\ ^ Ui + Si{vi{4>i) + \\i\). 

On the other hand, this line touches Ni{qs^-j(l)l''^-') only at the point {si,Vi{as^4>l^y), 
so that Ri{qsi-j(f>^'~''){y) is a constant and TrunCi(tp) does not divide qsi-j4>^^"'' ■ 



Ui + Si{vi{(i)i) + \Xi\) 



Therefore, Proposition 12 .11 shows that 

v{qs^-,{e,)cj,i{e,y^-^) = 

ei • • -Ci-i 
By the Theorem of the polygon |HN08[ Thm.3.1], 

Vi[(i)i) + |Ai| 



v{(t^i{e^)) 



Ui + j{vt{(t>i) + \h\) 



H, 



ei 



•ei-i 



(9) 

so that 

vfc.-j(6lp)) -- 

By©, ^;p(fe-,(0))=e(p/p)ff.J 
If e(p/p) = 1, then TTp := p is a pseudo-generator of p. If e(p/p) > 1, we can 
always find a pseudo-generator of p by computing a suitable product of quotients 
in the lists Quotj, divided by a suitable power of p. 



D 
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Corollary 3.3. 

(1) Let j'l, . . . ,jr be subindices satisfying, < ji < ei, for all 1 < i < r. Then, 
the following element belongs to "Zk ■' 

(2) If e{p/p) > 1, there is a unique family ji,...,jr as above, for which 
Wp(7rji,....jv) — 1. This family may be recursively computed as follows: 

jr = h;:^ (mod Cr), 
resr :~ {jrhr — l)/er, 

jr-1 ^ -h^\{Ur + jrVr{4>r) + reSr) (mod Er-l), 
reSr-1 := {jr-lK-l + "r + jrVr{4'r) +reSr)/er_i, 

ji = -h^^{u2 +J2V2i4'2) +res2) (mod ei). 

Proof. Item 1 is an imm.ediate consequence of item 1 of Proposition 13. 21 
Also, by Proposition 13.21 

Thus, item 2 states that there is a unique family ji, . . . ,jr such that 

1 
eip/p) 
Since e{p/p) = ei • • • e,. and \Xi\ — hi/ci, this is equivalent to: 

jlhi + U2 + i2V2{4>2) , 



-ffi.ji H 1- Hr,j,. = -j-rnx (™°^ ^)- 



Ml + jlt^l(0l) + - 



ei 



(mod Z). 



Clearly this congruence has a unique solution ji, . . . , jr satisfying < ji < ei, for 
all 1 < i < r, and this solution may be recursively obtained by the procedure 
described in item 2. D 

The only property of tp that we used in in Corollarv 13.31 is: ei ■ • • Cr = e{p/p). 
Thus, we don't need to use all levels of tp to compute a pseudo-generator of p; in 
practice we take r to be the minimum level such that ei • • • e^ = e(jp/p). 

3.2. Generators of the prime ideals. Let p be a prime number, and *p the set 

of prime ideals of K lying over p. Once we have pseudo-generators TTp G p of all 
p G *P, in order to find generators we need only to compute a family of integral 
elements, {bp £ Z^f }pgfp, satisfying: 

(10) «p(6p) = 0, Vpe'P, fq(6p)>l, Vq,pe*P, q^p. 

Then, for each p e *P, the integral element: 

ap := 6p7rp + ^ b^ e Zk 

clearly satisfies: Wp(Q:p) = 1, Wq(ap) = 0, for all q 9^ p. Therefore, p is the ideal 
generated by p and ap. The rest of this section is devoted to the construction of 
these multipliers {bp}. 
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Let (tp)pg(p be the parameterization of the set *p by a family of /-complete types 
obtained by an application of Montes algorithm. As usual, we suppose that each 
tp has been conveniently enlarged with an (rp + l)-th level, as indicated in section 
11.41 From now on we provide the invariants of tp with a subscript p to distinguish 
the prime ideal they belong to: rp, 0i^p, m^^p, A^^p, etc. 

The integral elements 6p will be constructed as suitable products of (/)-polynomials 
divided by suitable powers of p. The crucial ingredient is Proposition l3.8[ that com- 
putes Vp{(j)i,q{9)) for all p ^ q in *p, and all 1 < i < rq + 1. 

Definition 3.4. For any pair p, q G *P, we define the index of coincidence between 
the types tp and tq as: 



*(tp, tqj — 



0, iflpO-.p 7^ V'O.q, 

min{j e Z>o | (0j,p, Aj,p, ^j,p) J^ (0j,q, Aj,q, TZ-j^q)} , ifijQ.p = V'O.q- 

Alternatively, i(tp,tq) is the least subindex j for which TrunCj(tp) ^ TrunCj(tq). 
Remark 3.5. By definition, 

i>t,p = V'i.q, Vi < i(tp,tq). 



(^j.p 



■^i.p — Ai_q, 



Hence, by the definition of the p-adic valuations Vi^p, Vi,q, and by [HN08[ Thni. 
2.11], we get: 

Vi,p^Vi.q, m.j,p = mi,q, Wi,p(0i,p) = Wi,q(0j,q), Vi<z(tp,tq). 

Lemma 3.6. //p,qG*P, and p ^ q, </ien i(tp, tq) < minjrp, rq}. 



Proof. Suppose rp < r^ and i(tp,tq 



1. Then, TrunCrp(tp) = TrunCrp(tq 



Since the type Trunc,.p (tp) is /-complete, it singles out a unique p-adic irreducible 
factor of f{x) (item 2 of Definition II. ip . Hence, Trunc^p (tq) is also /-complete and 
it singles out the same p-adic irreducible factor of f{x). This implies that p = q. □ 



By dll) and 
(11) 



for all 1 < i < 



vp{4>^^p{e)) =e(p/p) 



1 we have 

Wi,p('/'j,p) + |Ai,p| 



ei,p 



^i,p 



In order to compute the values of i'p(0i,q(6')), for p 7^ q, we need still another 
definition. 



Definition 3.7. Let p,q £ 'P, q 7^ p, and j = i(tp,tq). Let Sp = #RefinementSj_p, 
and consider the list Refp obtained by extending the list RefinementSj^p by adding 
the pair [4'j,pi^j,p] o.t the last position: 



Ref„ 






{Sf) ^(Sp) 



'^J.p 



J,p 



(sp + 1) ;^(sp + l) 



'^3,P 



J,P 



i't'j.P ' Aj.p] 



Let Ref q be the analogous list for the prime ideal q . 

We define the greatest common 0-polynomial of the pair (tp,tq) to be the more 
advanced common (p-polynomial in the two lists Refp, Refq. We denote it by: 

</.(p,q):=0(5=</.(3, 

for the maximum index k such that <f>\ „ = 6) I . 

We define the hidden slopes of the pair (tp,tq) to be: Ap := A^ „, Aq := A^ q . 
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Remarks. 

(1) By the concrete way the processes of branching, enlarging and/or refining 
were defined, this polynomial (/)(p, q) always exists. In fact, let us show that we must 
have (JK I — 4>^ „. Since TrunCj_i(tp) — TrunCj_i(tq), this type had some original 
representative (say) (j)j. By considering NT „.(/) and the irreducible factors of 
all residual polynomials of all sides, we had different branches (A, ip) to analyze; if 
there was only one branch, the algorithm necessarily performed a refinement step, 
because otherwise we would have j(tp,tq) > j. After eventually a finite number of 
these unibranch refinement steps (that were not stored in the list RefinementSj ) , we 
considered some representative, let us call it 0^ again, leading to several branches. 
One of these branches led later to the type tp and one of them (maybe still the same) 
to the type tq. If the p-branch experimented refinement, the list RefinementSj.p had 
^ — 't'j as its initial ^-polynomial; if the p-branch was /-complete or had to be 



j,p 



enlarged, then the list RefinementSj^p remained empty and we had 
any case, 0, is the first (/)-polynomial of the list Refp. 



"hP 



U.J. 



In 



(2) All 0jp, 0^- q are representatives of tj_i :— TrunCj_i(tp) = TrunCj_i(tq); 
in particular, all these polynomials have degree rrij. With the obvious meaning for 
i/";- p, we have necessarily: 






x(fe) ,/,(fe) 



^ 






\(fe) ,/,(fe) 



C-A 



J\p' 



^],p 



Ai) ^(f) .M 



for all 1 < ^ < A;. Thus, 4>{p,c{) is the first representative of tj_i for which the 



branches of tp and tq are different. 



Ap 



(3) Caution: we may have Refp — Refq. In this case (/)(p, q) 



(^j.p 



/^j^q 



and 



A 



j.p 



A, 



Aq; the branches of tp and tq are distinguished by ipj^p ^ ipj.q- 



Proposition 3.8. Let p, q G *P, p ^ q, and j — i(tp,tq). Let (/)(p, q) be the greatest 
common (f>-polynom,ial of the pair (tp,tq) and Ap, Aq the hidden slopes. For any 
1 < i < rq + 1. 



vp{4>^,,m 

e(p/p) 



0, 



Viicpi) 


+ 


A.I 


ei •• 


ei 

+ 





ei • • -ej^i 
«,(0,)H-min{|Ag|,|AP|} 
ei • • • ej-i 
^,(0,)-hmin{|A^|,|AP|} 



H,q 'JjWj 



ei ■ ■ -e,-! 



«/ -i < j, 

ifi^j and ^j,q =(/)(p,q), 

ifi^j and(j>j^^ ^ 0(p,q), 

ifi>j>0. 



In these formulas we omit the subscripts p, q when the invariants of the two 
types coincide (cf. Remark \3.5\) . 

Proof. The case j = was seen in ([8|). The cases i < j and i = j, (j)j,q = 0(p, q), 
are a consequence of pTJl . 
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Suppose z > j > and 



'^3,P 



we have then, (/)(p,q) = 



'-'j.p 



and 



^p = ^j.p7 -^q = ^j,c[- We compute Vp{(f>i_q{9)) by applying Proposition 12.11 to the 
polynomial g{x) — 4>i^^(x) and the type tp. Since Vj,p — Vj,q and 0j^p — 0j^q, we 
have iVj.p(0i_q) = A'j_q(0i^q). On the other hand, we saw in the proof of Lemma FOl 
that iVj,q((/)i_q) is one-sided of slope Xj^q. Figure 3 shows the three possibilities for 
the line Lx^ ^ of slope Xjp that first touches iVj_p(0j,q) from below. 



h'\ 


, 




\\^j,q 


^ 


^^..X, 


','iJ 







mi^d/i^j 



H'-' 


. 


\ 


\^j,q 


2.q^ 


■'■■■■^^,,f 





mi^q/rrij 



1 

■(d)- 1 


:. \^j,q 

""^ T, ^ 




•••■■.:r.^j.p 





mi,(|/mj 



^i,p < '^3,^ 



^3:P — ^j:q 
Figure 3 



A glance at Figure 3 shows that 



H = Vj{4>i,q) + 



{\^],vl l^i,ql} = — ^ (t^i(0j) +min{|Aj-p|, |Aj,q|}), 



the last equality by Lemma [1.21 Now, if A^-^p 7^ Aj^q, the line L\^ ^ touches the 
polygon only at one point, and the residual polynomial Rj,p{4'i.c\)iy) is a constant. 
If Xj^p = Aj,q, then Lx^^ contains A'j\p((/'j,q) and i?j^p((/)j,q) = i?j,q(^i,q) is a power 
oitjjj^q, up to a multiplicative constant. In this case, necessarily ipj^q ^ V'j.p) by the 
definition of j = J(tp, tq). Therefore, TrunCj(tp) never divides 4'i,qi ^^d Proposition 
12.11 shows that v{(j)i^q{6p)) = H/{ei ■ ■ ■ ej_i). By ^, we get the desired expression 
for vp{(t)i^q{6)). 

Suppose now i > j > 0, (j)j,p ^ 4>j,qi or i — j, 0j_q 7^ 0(p, q). Consider a new type 
tp, constructed as follows: if 0j,p = (j){p, q), we take tp ~ tp, and if (/)j_p ^ 0(p, q), 
we take 



tp = ((/'i;Ai,(?!)2r 



,_i;Aj_i,(/.(p,q);A2, 



P' 



-^J,P ~ '^pj0J + l,p;-^j"+l,pi ■ • • ,0rp+l,p; Arp+l,p,'0rp+l,p 



By |GMN08[ Cor. 3. 6], tp is a type, and it is also /-complete. If (/)j,p ^ 0(p, q) then 
tp is not optimal because deg0(p,q) = deg^j^p, but optimality is not necessary to 
apply Proposition 12. II 

We consider an analogous construction for tq. The new types satisfy i(tp, tq) = j 
and they have the same j-th (/>-polynomial; finally if i = j, the polynomial (j)i,q is 
the (_7-|-l)-th 0-polynomial of tq. Therefore, the computation oivp{(f>i,q) is deduced 
by the same arguments as above. D 

We are ready to construct the family {&p}pe<p. Consider the following equivalence 
relation in the set *P: 

p - q ■^=^ ^o.p == -^o.q, 

and denote by [p] the class of any p G ^. 

For each class [p], let ^i,[p](x) G Z[x] be the first (/)-polynomial in any list Refq 
for some q G [p] ; we saw in the first remark following Definition 13.71 that all lists 
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Rcf q , for q G [p] , have the same initial (/)-polynomial. Now, for each q G [p], denote 
by A]* . the first sfope in the hst Refq. In other words, for any q G [p], we have 

I {4>i,q, ^i.q), if Refinementsi^q is empty, 
I [4>i q^^i a] ^ i^ Refinementsi^q is not empty. 
By © and ®, 

(12) -,(<^ifpl(e)) = | "' '''^^^^' 

^ ' '^"^'^^'^^ " \ e(q/p)|A?,J, ifqG[p]. 

Consider now, for each class [p]: 

Fix a prime ideal p G *p. If #[p] = 1, the element bp = (i?[p](6'))2 satisfies (fTO)) 
already. Suppose now #[p] > 1. For all [ G [p], I ^ p, let (f>i — (j)ri+i,i be the Montes 
approximation to fi{x) contained in ti; consider the least positive numerator and 
denominator of the rational number Vp{(j)[{9))/e{p/p) (which has been computed in 
Proposition [ 



Vp{(t)i{9)) _ rn 



, gcd(n[,d[) = 1. 



e(p/p) di ' 

We look for an integral element of the form: 



(13) hp ^ , 

where the exponents N, m are given by 

f iVe(q/p) + 2 



N = 2, '^Ii '7^ 

I6[P],I#P 

Take q ^ [p]. By dH) and ([12]), we have 



qeqJ.q^M I e(q/p)|AO J 



n '^'(^)'' = 0' «q(^[P](^)) - «q('^i,M(^)) - e(q/p)|AO,J. 

je[p]-i5^p / 

Hence, Vq{p^bp) — me(q/p)|A'j' ^j > Ne{q/p) + 2, so that Uq(6p) > 1, as desired. 
For the prime p itself, we have Vp{B!pi{9)) — and, by construction, 

vpibp)=[ E d,Vp{M9))]-Ne{p/p) = 0. 
\ie[pl.i#p / 

Finally, for a prime ( G [p], [ ^ p, we have vi{B[p-^{6)) — and 

vi{bp) = Vi + V2~Neil/p), Vi:= ^ vii^vi9f^'), V2 -.^ Vi{MOf')- 

i'e[p],i'#p,[ 

Lemma 13.61 and Proposition 13.81 show that Vi (as all invariants we used so far) 
depends only on the numerical invariants of the types t(,t[', of level 1 < i < 
*(tt,t(') < min{r[,r['}, and not on the quality of the Montes approximations (pv. 
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On the other hand, V2 depends on the choice of (j)i as a Montes approximation of 
Mx); by m- 

ei,f • ■er^,l 
Hence, for ah [ G [p], I 7^ p, we improve the Montes approximation 0( till we get 

2 + Ne{\/p) - V, 

hri + lA > 1 Vri + l,l[<Pl). 

"I 

This ensures that vi{bp) > 1, as desired. 

3.3. Two-element representation of a fractional ideal. Any fractional ideal 
a of ii" admits a two-element representation: a = (^, a), where a ^ K and £ = 
£{a) e Q is the least positive rational number contained in a. It is straightforward 
to obtain such a representation from the two-element representation of the prime 
ideals obtained in the last section. For the sake of completeness we briefly describe 
the routine. For each prime ideal p of K, we have computed an integral element 
ctp € Zk such that 

wp(ap) = 1, 'yq(ap) = 0, Vq|p, q T^p. 

These elements are of the form: ap — p''''hp{6), where i^p G Z and hp{x) G Zlx] is 
a primitive polynomial. Generically, J^p < 0, except for the special case p = pIjk, 
where i^p = 1, hp{x) — 1. Let us write 

NK/Q{hp{e))^p>'''Np, withpfiVp. 

Suppose first that o — IlpinP"'' ^^^ support only in prime ideals dividing p. We 
take then: 

" = n"p' n ^P°''' ^ '■= maxj- 

PIp ap<0 '^'^ *-' 

One checks easily that £{a) = p^ and a = {p^ , a). 

In the general case, we write o = ripep '^p^ where P is a finite set of prime 
numbers and a^ is divided only by prime ideals lying over p. For each a^ we find 
a two-element representation Op — {p^^ ,ap)\ then, the two-element representation 
of a is: 

^^p"^El n '^"'^M«p 

^pev pev \qev,q^p 

Note that the second generator a constructed in this way satisfies: Vp{a) — Vp{a), 
for all p with Vp{a) ^ 0, which is slightly stronger than the condition o = {£,a). 

4. Residue classes and Chinese remainder theorem 

In this section we show how to compute residue classes modulo prime ideals, and 
we design a Chinese remainder theorem routine. As in the previous sections, this 
will be done without constructing (a basis of) the maximal order of K and without 
the necessity to invert elements in the number field. Only some inversions in the 
finite residue fields are required. 

For any prime ideal p oi K we keep the notations for tp, fp{x), (j)p{x), Fp, 9p, 
Kp , "Lk^ , as introduced in section [2j 



. e(p/p) 
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4.1. Residue classes modulo a prime ideal. Let p be a prime ideal of K, 
corresponding to an /-complete type tp with an added (r + l)-th level, as indicated 
in section 11.41 

The finite field Fp := F^+i may be considered as a computational representation 
of the residue field "Lk /p- In fact, fix the topological embedding, ip-. K ^^ K^, de- 
termined by sending 6 to 9p, and consider the reduction modulo p map obtained by 
composition of the embedding Zk ^^ ^A'p with the local reduction map constructed 
in dlD: 

Ircdp 

redn : Zk '^ Zk„ — > ¥„ 



ip: ^K ^^ ^A'p — ^ J^ p ■ 



The commutative diagram: 



; ; ^ II 

Zk/p -^ ^kJpI-k, ^ Fp 

shows that our reduction map redp coincides with the canonical reduction map, 
Zfc — > Zif/p, up to certain isomorphism Zk/P ^^ Fp. 

The problem has now a computational perspective; we want to find a routine 
that computes redp (a) G Fp for any given integral element a G Zk- To this end, 
it is sufficient to have a routine that computes lredp(Q;) G Fp, for any p-integral 
a G K. Let us show that this latter routine may be based on item 2 of Proposition 
[lTT]and Lemma [Lil 

Any a G Zk can be written in a unique way as: 

ag{0) 
cy ^ — 

b p^ ' 

where a, b are positive coprime integers not divisible by p and g{x) G Z[x] is a 
primitive polynomial. Clearly, 

redp (a) = Ircdp (tp (a)) = lredp(a/6) lredp(5(6'p)/p^), 

and Iredp (a/6) G Fp is the element in the prime field determined by the quotient of 
the classes modulo p of a and b. Thus, we need only to compute Iredp (g(0p)/p^). 

If A^ = 0, then lredp(5(6'p)) — iedp{g{6)) is just the class of g{x) modulo the 
ideal (p, (/)i.p(x)). In other words, ii g{x) G Fo[x] is the polynomial obtained by 
reduction of the coefficients oi g{x) modulo p, then Tedp{g{9)) = g{zo) G Fi_p C Fp. 

If TV > 0, we look for the first index, 1 < i < r -I- 1, for which the truncation 
TrunCi(t) does not divide g{x). In the paragraph following Proposition 12.11 we 
showed that this will always occur, eventually (for i = r + 1) after improving the 
Montes approximation (/)p = 4>r+i- By Proposition 12.11 there is a computable point 
is,u) eN,{g) such that u(g(6'p)) = {sh, + ua) / {a ■ ■ ■ e^) = w($,(6'p)"7r,(6lp)"), and 



9(9p} 
$z(ep)^vr,(0p)'^ 



Iredp ( ^ ,/:Z\, ,„ ] = i?.(5)(^.) ^ 0. 



Now, if {shi + uei)/{ei ■ ■ ■ Ci) > N , we have Iredp (.g(6'p)/p^) = 0; on the other hand, 
if (shi + uei)/{ei ■ ■ ■ e,;) — N, we have 

(14) Iredp (^) = Rdg){z.) ■ hedp (^ill^^ZjlM!) ^ i?^(,)(,,) ,*.... ,*., 
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where p^^$i(6'p)'*7ri(0p)" = 7*^ • • • 7'% and the vector (ti, . . . ,ti) G Z' can be found 
by the procedure of Lemma 11.41 apphed to the input vector 

log (p-^$,(x)^7r,(x)") = {-N, 0, . . . , 0) + s log$, + ulog^,. 

Since log$i, logTr^ have been stored as secondary invariants of tp, we get in this 
way a really fast computation of lredp(p^^g(0p)). 

This ends the computation of the reduction modulo p map redp . 

4.2. Chinese remainder theorem. It is straightforward to design a Chinese re- 
mainders routine once the following problem is solved. 

Problem. Let p be a prime number, *P the set of prim.e ideals of K lying above p, 
and (ap)pg(p a family of non-negative integers. Find a family (cp)pgq3 of integral 
elements Cp G Iik such that, 

Cp = 1 (mod p"^), Cp = (mod q""), Vq G <P, q^p, 

for all p G *p. 

There is an easy solution to this problem: take the element 6p G Z/f satis- 
fying ([To]) , constructed in section [3?2| and consider Cp = [b^Y ('p^^', where qp — 
N^/q(P) = #Fp, and t is sufficiently large. However, the element Cp constructed 
in this way is not useful for practical purposes because it may have a huge norm 
and a huge height (very large numerators or denominators of the coefficients of its 
standard representation as a polynomial in 6), if q^ or t are large. Instead, we shall 
refine the construction of the element 6p to get a small size solution Cp to the above 
problem. 

First we deal with the particular case Op = 1. The idea is to get an element 
6p G l^K satisfying 

Wp(6p) = 0, 'yq(5p) >aq, VqG*P, qT^p, 

and then find /3 ^ K such that Vf,{f3) = 0, Cp := bp/3 is integral, and Cp = 1 (mod p). 
Since we know the two-element representation q = {p,aq) of the prime ideals, we 
could take 6p — Yiaev q=^p('^ti)"'' ■ Again, this might lead to a 5p with large size, 
so that a direct construction of bp is preferable in order to keep its size as small as 
possible. 

Thus, we consider an element bp G Zk as in (113^ : 

, iB[p]ix)rUi^[p},i^pM^r' 

h = ^^v ' 

with N = E[e[p], MP div{(l)i{0p)) = Eie[p], Mp '^'- ^^ construction, Vp{bp) = 0. 

Let i = maxqgrp q-^p{i(tp,tq)}; Lemma [3.61 shows that i < rp. From now on, we 
use only invariants of the type tp and we drop the subindex p in the notation. Let 



0, if i = 0, 



M 



ei ■ • • Ci 



if i > 0. 



Arguing as in section [3.21 we can take m sufficiently large, and each (/"[(a;) suffi- 
ciently close to the p-adic irreducible factor fi{x), so that 

(15) vq{bp) > aq + A/e(q/p), Vq G ^3, q ^ p. 
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while keeping the denominator p^ and the condition Wp(6p) = 0. In particular, 
hp belongs to Iik- The idea is to multiply bp by some element in K that conve- 
niently modifies its residue class modulo p. We split this task into two parts, that 
may be considered as a kind of respective inversion modulo p of (_Brpi(0))™ and 

Let h{x) = (B[p](x))"'. Then, C := redp(/i(6l)) G Fi C Fp is just the class 
of hix) modulo the ideal {p,(j)i{x)). We invert C, in Fi and represent the inverse 
C,^^ = P{zo), as a polynomial in zq of degree less than /o, with coefRcients in the 
prime field Fq. Take /3o = P{d), where P{x) G Z[a:] is an arbitrary lift of P{x); 
clearly, /i(0)/3o = 1 (mod p). 

Let now gix) = P'"" Ilielp]A^pMxY' ■ K [p] = {p}, then i = 0, N = M = 0, 
g{x) = 1 and we are done. Suppose [p] 2 {p}: so that 1 < i < rp. By the definition 



of the index of coincidence, we have TrunCi(t 
of the product [HN08I Thm.2.26], TrunCi(tp 
last section. 



p)\M^) 



\9{^)- 



e:=lredp(p-^\g(0p)) = i?,(g)(z,)- 



for alH 7^ p; by the Theorem 
Therefore, as we saw in the 



gFh 



for some easily computable sequence of integers {ti, . . . ,ti). Let V — ei ■ ■ ■ e^M; by 
[HN08[ Cor.3.2], v (7r,+i(6'p)^) = M. Compute a vector {t[,. . . ,i^) G Z* such that 



p-^7r,+i(x)^ 



7i 



■ 7j', as indicated in Lemma ri.4[ and take 



e--^ 



t[ 



t' 



eF,+i 



Let if{y) G Fj[y] be the unique polynomial of degree less than fi, such that ^'(zi) — 
e 

Zi 



J.vie,^^^yx_ Let ^ = ordy(^(zj). Clearly, 

y > Vi+i{(l)i+i) = e.JiV.i+i{ 



0, 



the last equality by [HN08[ Thm.2.11]. Therefore, we can apply the constructive 
method described in [HNOSl Prop. 2. 10] to compute a polynomial P{x) G Zlx] sa- 
tisfying the following properties: 

degP(x)<m,+i, v^+i{P) = V, y''R,{P)iy)^^iy). 

A look at the proof of [HNOSi Prop. 2. 10] shows that Ni{P) is one sided of slope Xi 
and its end points have abscissa CiV and ej deg f (cf. Figure 4). 




Figure 4 
Claim: lredp(p-^P(6lp)) = C"^ 

In fact, since degP(a;) < rrii+i and Vi+i{P) = V, the Newton polygon Ni+i{P) 
is the single point (0,V^); in particular, TrunCi+i(tp) does not divide P{x) and 
Proposition 12.11 shows that v{P{9p)) — M and 



lredp(P(0p)A,+i(^p 



P.+i(P)(z.+i)^0. 
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Actually, Ri+i{P){y) has degree and it represents a constant in F^+i that we 
denote simply by Ri+i{P). By (fM]) . 

Iredp {p-''P{e,)) = B.,+i{P) lredp(p-*^^,+i(0p)^) = i?.+i(P)C'. 

By the very definition of the residual polynomial (HN081 Defs.2.20+2.21], we have 

and the Claim is proven. 
Finally, consider 



c 



= b,MP'''P{9)) 



p — "p 

The condition (fT5|) ensures that Cp belongs to "Ek (although p^^^ P{9) might not 
be integral) and satisfies i'q(cp) > a^, for all q G *P, q 7^ p. By construction, 
redp(cp) = Iredp (ip(cp)) = 1. 

It remains to solve our Problem when Op > 1. In this case, we find c G Z^ such 
that 

c= 1 (modp), ceeO (mod q°'i), Vq e«P, q =^ p, 

and we take Cp = (c — 1)'" + 1, where m is the least odd integer that is greater than 
or equal to ap/up(c — 1). 

5. p-INTEGRAL BASES 

Let p be a prime number and let Vp be the prime field of characteristic p. 
A p-integral basis of K is a family of n Z-linearly independent integral elements 
ai, . . . , a„ G Zif such that 

p\ (Zk: {ai,...,an)^) , 

or equivalently, such that the family ai ® 1, . . . , a„ 1 is Fp-linearly independent 
in the Fp-algebra Zk ®i Fp. 

If the discriminant disc(/) of f{x) may be factorized, the computation of an 
integral basis of K (a Z-basis of Z/f ) is based on the computation of p-integral 
bases for the different primes p that divide disc(/). 

Anyhow, even when disc(/) may not be factorized, the computation of a p- 
integral basis of K for a given prime p is an interesting task on its own. In this 
section we show how to carry out this task from the data captured by the Okutsu- 
Montes representations of the prime ideals p lying over p. For any such prime ideal 
we keep the notations for fp{x), 4'p{x), Fp, 9p, Kp, I^Kp, as introduced in section [2] 

5.1. Local exponent of a prime ideal. Let *p be the set of prime ideals lying 
over p. For any p G *P we fix the topological embedding K ^-)- Kp determined by 
sending 6* to 6'p. 

Definition 5.1. We define the local exponent of p E ^ to be the least positive 
integer exp(p) such that 

Note that exp(p) is an invariant of the irreducible polynomial fp{x) G Zp[a;], but it 
is not an intrinsic invariant of p. 
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The computation of exp(p) is easily derived from the resuhs of I GMN09b] . 
Let p = [p;ipi, . . . ,(j)r, (j)p] be the Okutsu-Montes representation of p, as indicated 
in section O By |GMN09b[ Lem.4.5]: 

fp{x) = (t>p{x) (mod m^'^T), where u ^ Vp + (hr+i / e{p / p)) , 



and Vp is the rational number 



hi 



ei 6162 ei • • -Cr 



Caution: this number is not always an invariant of fp{x). If e^fr = 1, the Okutsu 
depth of fp{x) \s R~ r — 1 (cf. ([3])), and hr depends on the choice of 4>r- 

Denote 4>o{x) — x, mg — 1, m^+i = rip := e{p/p)f{p/p). Any integer < to < 
rip can be written in a unique way as: 



rrii+i 



^fljirij, < aj < 



i=0 

Then, g„i (x) :— nl=o 't'i (^)"' ^^ ^ divisor polynomial of degree m of /p (x) [GMNOQbl 
Thm.2.15]. Also, if i/„i := lv{gmiO))\, the family 

is a Zp-basis of 'ZiKp [Oku821 I,Thm.l]. Since the numerators gm{x) have strictly 
increasing degree and J^i < • • • < J^np-i, it is clear that exp(p) = Vn^-i- 
On the other hand, since mi^i/rrii — Cifi, we have 

r 

Now, along the proof of [GMNOQbl Lem.4.5], it was proven that: 

(17) ±{e.h ~ IMMO,)) - vi^^^M) - -P - ^ = "'ll^/lf -P- 

Thus, if we combine (J17p with the explicit formula: 

r 
i=l 

given in jHN08| Prop. 2. 15], we get the following explicit computation of exp(p) in 
terms of the invariants of the Okutsu-Montes representation of p. 

Theorem 5.2. For all p e *p, we have exp(p) = [/ipj, where 

Vr+l{(t)r+l) V^. „ J- l\ ^« 

Note that \Xp is an Okutsu invariant, because it depends only on ei,fi,hi, for 
1 < I < i?, where R is the Okutsu depth of /p(a;). If i? = r — 1, then e^/r = 1 and 
the summand of /ip corresponding to r = r vanishes. 
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5.2. Computation of a p-integral basis. Along the computation of the prime 
ideal decomposition of the ideal 'pLk (by a single call to Montes algorithm) we can 
easily store the local exponents exp(p), and the numerators gmi^x) and denominators 
p"^™ of all Zp-bases of Z^p given in (fT6|) . for all p e *p. 

It is well-known how to derive a p-integral basis from all these local bases. Let 
us briefly describe a concrete procedure to do this, taken from IOre25j . 

We apply the method described at the end of section 13.21 to compute multipliers 
{6p}peq5 satisfying: 

(18) «p(6p) = 0, «q(6p)>(exp(p)-|-l)e(q/p), Vqe'P, q^p. 

Consider the family obtained by multiplying each local basis (J16p by its correspond- 
ing multiplier: 



^P 






Then, B :— Upeoj '^p ^^ ^ p-integral basis of K. 

In fact, although the elements gra(,G)lv^'^ £^re not (globally) integral, (IT51) shows 
that the products Q;„i,p := 6p (<?m(^)/p'''") belong all to TLk and satisfy 

^'q(am,p) > e(q/p), Vq G <P, q =^ p. 

It is easy to deduce from this fact that the family of all oim^p determines an Fp- 
linearly independent family of the algebra Ik ®% Fp. 



6. Some examples 

We have implemented the algorithms described above in a package for Magma. 
The arithmetic of number fields in any algebraic manipulator has to face two prob- 
lems: the factorization of the discriminant and the memory requirements for large 
degrees. Our package allows the user to skip the first problem, while it uses very 
little memory, expanding Magma's capabilities by far. The package which can be 
downloaded from its web page (http: /ma4-upc.edu/^guardia/+Id.eals.html), 
is described in detail in the accompanying paper [GMNlOj. We include here a 
few examples which exhibit the power of the package in different situations. More 
exhaustive tests of Montes algorithm have been presented in |GMN08| , |GMN09a] . 

The computations in these examples have been done with Magma v2.15-ll in a 
Linux server, with two Intel Quad Core processors, running at 3.0 Ghz, with 32Gb 
of RAM memory. 

6.1. Large degree. Consider the number field K — Q{9) given by a root 6* G Q of 
the polynomial f{x) = x^"°" + 2^°a::^° + 2^°. The factorization of the discriminant 
of/ is 

Disc(/) = 25394"3^"52"°°127^°313^°743^°4886229527^V", 

withp = 337572698551220494882323528404563236947916489629537. The large de- 
gree of / makes impossible to work in this number field using the standard func- 
tions of Magma, even after factorizing the discriminant, since the computation of 
the integral basis is necessary for these functions. But our algorithms avoid this 
computation, so that we can work with ideals in K. For instance, in the table 
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below we show the local index of the primes dividing Disc(/) and the time taken 
to decompose them in K . 



Ideal 


Index 


Time 


2Zk 


26235 


0.36s 


3Zk 





0.61s 


5Za- 


20 


0.63s 


127Za' 





1.29s 


313Za 





3.69s 


743Za 





6.47s 


4886229527Za' 





6.96s 


pZk 





60s 



Thus, we need less than 90 seconds to see that the discriminant of K is 

Disc(ii:) = 2i4™3^°5i^^"l27'^°313^°743^°48862295275V°- 

The running times in the table show clearly that the cost of the factorizations 
increases mainly because of the size of the numbers involved, and that the index 
has not a serious impact on them. The largest type appearing in these computations 
has order 3, and it appears along the factorization of the ideal 2Za:, which is 

2£k = Pi (Pij p4 (p4J p20, 

where p^ stands for a prime ideal with residual degree / and ramification index e. 
While we cannot expect to factor the ideals / = {6^ + 50)Zk, J = {6 + 10)Za: in a 
reasonable time, it takes 0.03 seconds to compute the factorization of its sum; 

I + J ^ pl{p[)Y^{p',)Y2o- 

The decomposition of 5 in the maximal order Za is 

^V — m5('«' ^5„20^„/ ^20„25„25„25('„/ n25 
5Za -p2(P 2) P2 IP 2) P2 P4 P15IP15J ■ 

With the residue map computation explained in subsection 14. 1[ we may check very 
quickly that 

6lEEC(modp4), 
where Za/P4 — IFsiC]) with ("* + 2(^ + 3 = 0. The Chinese remainder algorithm 
works also very fast in this number field. 

6.2. Small degree, large coefficients. The space of modular forms of level 1 
and degree 76 has dimension 6. The newforms in this space are defined over the 
number field K = Q{6), where 6* £ Q is a root of the polynomial: 

f{x) = a;6 + 57080822040x5 - 198007918566571424544768x4 
-11405115067164354385292006554337280x3 
+9757628454131691442128845013041495838774263808x2 
+290013995562379500498435975003716024800114593761580810240X 

-92217203874207784163935379997152082331434364841943058919508374716416. 

The discriminant of /(x) is 

Disc(/) - 2264372516716112^32^7419243259 . 1932. 

293 • 3919872475942721376796815212m, 

where M is a composite integer of 135 decimal figures which we have not been able 
to factorize. J. Rasmussen asked us ( [RaslOj ) for a test to check certain divisibility 
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conditions on the ring of integers of K ^ related to his work on congruences satisfied 
by the coefficients of certain modular forms. The time to find the decomposition of 
the primes in the set 

S := {2, 3, 5, 7, 11, 13, 17, 19, 43, 59, 193, 293, 391987, 4759427, 137679681521} 

is almost negligible, since it involves only types of order at most 1. The table below 
shows the local indices of these primes: 



Ideal 


Index 


2Za- 


132 


3Zk 


36 


5Za- 


8 


7Za- 


8 


IIZa 


1 


13Za 


1 


ITLk 


2 


19Zk 


1 


43Zk 


1 


59Zk 





193Zk 


1 


293Zk 





391987Za: 


1 


4759427Za 


1 


137679681521Zk 


1 



p2plp'iPiPi'. 



Hence the discriminant of K is Diac{K) — 59 • 293iV, where N is divisible by at 
least one of the prime factors of M, since M is not a square. The ideal prime 
decomposition of 3 in Zx is 

3Za' 

The algorithm explained in section [3] provides generators for all these ideals: 

P2 = 3Za+ 3-^^(46^+43116^+ 17170386|3+ 29006916|2+ 8201256*+ 2834352)Zif 
pi = 3Za+ 3-"(26i5 + 18156'4 + 5869806(3 ^ 7321596/2 + 65828761 + 1535274)Za 
20316'''+ 662796613+ 11236326i2- 

- 23076*4 + 9108726*3 + 8475846(2 

h 20916*4+ 7086966*3+ 6463806*2-1 

Applying the algorithm described in section 21 we can compute without much effort 
an element a € K satisfying 

l(mod P2), a = 6* (mod pi). 



P'l = 3Za- 
P'i' = 3Za- 
p'l" = 3Za 



3-^(26*54 

3-"(26*5- 

h3-"(26*^- 



- 10716306* + 295245)Zk 

f 3980346* +1121931Zx 

6342306*+ 1121931)Zif 



a 



a 



(mod(pi)^ 



3d (p'/)^) , a = 6* (mod (p;")4) 



We may take, for instance: 

a = 3^^(7860866*5 ^ 4459896*" + 1968576*3 + 11593536*2 + 6495396* + 354294). 

Following the algorithm for p-adic valuations introduced in section [21 we can check 
this result computing the valuations of the differences r — 6^ at the prime ideals 
dividing 3: 

t;p, (a - 1) = 1, «p, {a -6) = 7, Vp>^ {a - 6^) = 4, v^n [a - 6^) = 4, v^,n (a - 6^) = 4. 
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All these computations are almost im.mediate. Even the computation of an S- 
integral basis takes only 0.06 seconds. 

6.3. Medium degree. Consider the polynomials: 

00 = a; + 1, 01 = 0Q + 2, 

021 = 0? + 8, 022 = 0f + 4000? + 32, 

03 = 0^2 



^2 ^25602, / = 03021 +230. 



Let K = Q{0) be the number field of degree 20 determined by a root € Qof f. For 
the prime p = 2, the polynomial / has two complete types, with associated Okutsu 
frames [0i, 02i] and [0i, 022, 03], which give rise to the two prime ideals of Z/^ over 
2. The concrete decomposition is 2Zk — PiPI' where /(pi/2) = i, e{pi/2) — Ai. 
The discriminant of / is 

Disc(/) = 2268 . 32 ■ 19927 • 43691^ ■ 211039 • 6059454913- 

512920919154157817 • 25506978885046388417449- 
149169795543042282387542317948232968678925571739. 

In this example we may compare the performance of the standard Magma func- 
tions and that of our package, since Magma can determine the ring of integers of K. 
Once the factorization of Disc(/) is known, Magma takes 5.8 seconds to determine 
Z/f, and 0.08 seconds to find the decomposition of the prime 2 in Z^. Our package 
takes 0.3 seconds to see that Disc(iir) — 2~234 Disc(/), and during this computation 
already finds the decomposition of all the primes dividing the discriminant. Our 
program can also compute a 2-integral basis of K, which is already a global integral 
basis, in 0.02 seconds. 

7. Conclusions 

7.1. Challenges. We described routines to perform the basic tasks concerning 
fractional ideals of a number field, based on the Okutsu-Montes representations of 
the prime ideals [Mon99j , |HN08] . This avoids the factorization of the discriminant 
of a defining equation and the construction of the maximal order. These routines 
are very fast in practice, as long as one deals with fractional ideals whose norm 
may be factorized. 

A big challenge arises: is it possible to combine these techniques with some kind 
of LLL reduction to test if a fractional ideal is principal? 

Also, the generators of the prime ideals constructed in this paper have small 
height as vectors in Q" (the coefficients of its standard representation as a poly- 
nomial in 9). This may have some advantages, but in many applications it is 
preferable to have generators of small norm. A solution to the above mentioned 
challenge would probably lead to a procedure to find generators of small norm too. 

7.2. Comparison with the standard methods. Suppose the discriminant of 
the defining equation of the number field may be factorized. Most of the methods 
to compute a Z-basis of the maximal order are based on variants of the Round 2 
and Round 4 algorithms of Zassenhaus. The procedure of section [S] yields a much 
faster computation of an integral basis and the discriminant of the field. 

Once the maximal order is constructed, we can compare our routines for the 
manipulation of fractional ideals with the standard ones. The routines based on 
the Okutsu-Montes representations of the prime ideals are faster, mainly because 
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they avoid the usual hnear algebra techniques (computation of bases of the ideals, 
Hermite and Smith normal forms, etc.), which become slow if the degree of the 
number field grows. 



7.3. Curves over finite fields. The results of these paper are easily extendable 
to function fields. If C is a curve over a finite field, there is a natural identification 
of rational prime divisors of C with prime ideals of the integral closures of certain 
subrings of the function field lHess99j , |Hess02| . Montes algorithm may be applied 
as well to construct these prime ideals, and the routines of this paper lead to parallel 
routines to find the divisor of a function, or to construct a function with zeros and 
poles of a prescribed order, at a finite number of places. 

The results of section [5] may be used to efficiently compute bases of the above 
mentioned integral closures too. However, the big challenge of section 17.11 has its 
parallel in the geometric situation: we hope that the techniques of this paper may 
be used to find better routines to compute bases of the Riemann-Roch spaces and 
to deal with reduced divisors. This would open the door to operate in the group 
Pic (C) of rational points of the Jacobian of C, for curves with plane models of 
very large degree. 
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